Intrinsic Geometric Structure of c = —2 Quantum Gravity 

J.Ambj0rn a , K.N.Anagnostopoulos a , T.Ichihara b , L.Jensen a , N. Kawamoto , Y.Watabiki b and 
K.Yotsuji c 

a Niels Bohr Institute, Blegdamsvej 17 
2100 Copenhagen 0, Denmark 

b Tokyo Institute of Technology, O-okayama 
Meguro, Tokyo, Japan 

c Department of Physics, Hokkaido University 
Sapporo, Japan 

We couple c = —2 matter to 2-dimensional gravity within the framework of dynamical triangulations. We use a 
very fast algorithm, special to the c = — 2 case, in order to test scaling of correlation functions denned in terms of 
geodesic distance and we determine the fractal dimension du with high accuracy. We find du = 3.58(4), consistent 
with a prediction coming from the study of diffusion in the context of Liouville theory, and that the quantum 
space-time possesses the same fractal properties at all distance scales similarly to the case of pure gravity. 



1. The model 

Our starting point is the partition function of 
Euclidean 2d quantum gravity at hxed volume 
and discretize it by using dynamical triangula- 
tions with a fixed number of equilateral triangles 
N. The integral over metrics changes into a sum 
over triangulations TJv, and by coupling the quan- 
tum gravity to c Gaussian fields, the partition 
function is rewritten as 

Z N = J2^-(det'Cr N )- c / 2 (1) 

Here Ct n is the adjacency matrix of the (j> 3 
graph, which is dual to TJv, St n is a symmetry 
factor, and det'Cq- N denotes the determinant of 
Ct n with the zero-eigenvalue removed. This par- 
tition function also serves as the defining equa- 
tion for coupling to matter with c < 0. We re- 
strict ourselves to surfaces with the topology of 
a sphere, for reasons which will soon be clear. 
By selecting c = —2 it is possible to construct 
a recursive sampling algorithm for the numerical 
generation of an ensemble of surfaces ||. Two 
results from graph theory are crucial here. First, 
the determinant (appropriately defined) of Ct n 
is equal to the number of spanning trees in the 



4> 3 graph. Second, a graph can be embedded on a 
sphere if - and only if - it is planar. Thus restrict- 
ing ourselves to configurations with the topology 
of a sphere, we can decompose the sum over all 
triangulations into a sum over all combinations of 
4> 3 trees and rainbows - a rainbow being a set of 
non-crossing lines which connect the endpoints of 
the tree-graphs. Both the trees and the rainbows 
satisfy Schwinger-Dyson equations, and by solv- 
ing these equations, we can generate the correct 
ensemble of configurations using a recursive sam- 
pling algorithm. The time for generating a con- 
figuration grows linearly with its size, a rare case 
of a non-trivial statistical system with subexpo- 
nential time growth of statistically completely in- 
dependent configurations. This fact has enabled 
us to generate a substantial number of configura- 
tions with up to 8 million triangles. 

2. Definition of the fractal dimension 

The fractal dimension is defined from scaling of 
correlation functions defined in terms of geodesic 
distance r, therefore we have to define geodesic 
distance on the triangulated surfaces. This is 
done in two ways. Either as the number of links 
between vertices in the cf> 3 graph (the dual lat- 



tice) , or the number of links between the vertices 
of the triangles (the direct lattice). It is a non- 
trivial fact that both definitions are proportional 
to each other after taking the quantum average. 
In order to study the intrinsic structure, consider 
the average volume n n {R) of a spherical shell at 
a distance R. We can then define two fractal di- 
mensions dh and dn in the following way: 



n N {r)=N 1 - l ' dH F 1 {x), 



x = 



r + a 

J\fl/d H ' 



(2) 



where Fi(x) <~ a; dh_1 for x <~ 0. The existence of 
the scaling variable x at all length scales is a non- 
trivial fact to be tested by the simulations. The 
"shift" a is a finite size correction. If space-time 
looks the same at all scales then dh = da- Our 
simulations provide evidence, that this is indeed 
the case. For the correct value of dn and a, the 
distributions nnir), plotted as a function of x, 
will fall on top of each other if they actually do 
scale. That this is indeed the case, and that the 
finite-size correction a is necessary, can be seen 
in fig. 1 and fig. 2 were the dual lattice distance 
has been used and the configuration sizes range 
from 2k to 8M triangles. 




Figure 2. F^x) with shift (a = 4.50) 
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Figure 3. Short distance higher moments 
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Figure 1. Fi(x) without shift (a = 0) 



3. Higher moments 

We can also consider the higher moments 
(l n )r,N of the boundary length I. These are 
defined in terms of the loop-length distribution 
function p(r,l), which counts the number of dis- 
connected loops of the set of points at distance 
r from a given point, as (l n ) r% N = J dll n p(r,l). 
From Liouville theory we expect [I 2 ] = [A], but 
the simulations are consistent with the scaling 

(l n ) r , N = N^F n (x) ,n>2, (3) 

where F n (x) <~ x 2n for x <C 1. This implies that 
[; d «/ 2 ] = [A]. The small distance scaling can be 
seen clearly in fig. 3 where we show a plot of the 
logarithmic derivative of the moments n = 2, 3, 4. 
Trying to collapse the distributions of the mo- 



Table 1 

The fractal dimension determined from different methods. 







lattice 


method 


dh 


3.58(5) 


direct 


Short distance scaling 


d h 


3.50(4) 


dual 


Short distance scaling 


d H 


3.60(3) 


direct 


Collapsing the 2-point function 


d H 


3.55(3) 


dual 


Collapsing the 2-point function 


d H 


3.577(8) 


direct 


Intersection of R a ,N 


d H 


3.574(20) 


dual 


Intersection of R a ,N 


d H 


3.64(7) 


direct 


Collapsing of higher moments 



ments is consistent with eq.(Q). 

4. Determination of fractal Dimension 

Using different methods to collapse the distri- 
butions jy we can determine both dn and dh 
with great precision. The results using these 
methods are shown in table 1. R a ,N(dii) = 
l j N i+x/dH j2 r (r + a)n N (r) is the average radius 
of the universe JlJ and dn can be determined by 
tuning du and a so that the value of R a ,N{dH) 
is independent of N. As is shown in fig. 4 with 
the link distance, this works very well indeed. 
Two conclusions stand out. First, dh and dn 
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Figure 4. R a ,N on the direct lattice. 



are equal - i.e. the generated surfaces look the 
same at small as well as at large distances. This 
is a non-trivial observation which is not always 
satisfied, for example in multi-critical branched 
polymers KM- Then, there exists two different 



predictions for the value of dn- One is obtained 
by using the diffusion equation in Liouvillc the- 
ory || which predicts 



= 2 V25^_T V ^_ , = 

V25^+VT^ 

whereas redefining the distance in the context of 
transfer matrix theory[Q gives 

2 24 



dn 



l7(c)| 



1 



2 (5) 



c+y/(25-c)(l-c) 

where 7(c) is the string susceptibility. Compar- 
ing these with the results in table 1, it is apparent 
that the Liouvillc-prcdiction is strongly favored, 
whereas the prediction from transfer matrix the- 
ory is disproved. 
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